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Abstract 

Different kinds of instabilities (CDW, SDW, SS) in the ID Hub- 
bard model with pair-hopping interaction are investigated using an ap- 
proximate Bethe-Salpeter equation. The study is performed at any den- 
sity of electrons and for arbitrary values of the model parameters. In 
the absence of the on-site interaction, no transition occurs at half filling 
for any finite negative pair-hopping parameter, in agreement with re- 
cent results; our calculation suggests that the Penson-Kolb model (t, W) 
with | W/t |< 7r/ sinfcp behaves qualitatively like the Hubbard model 
(t, U). Phase diagrams of the Penson-Kolb-Hubbard model (t, U, W) at 
various densities are presented. 
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1 Introduction 



The studies of the one-dimensional models of correlated electronic systems 
can be a primary source to understand the occurence of the high temperature 
superconductivity in materials which physics is mainly two-dimensional. In 
these high-T c superconductors the "Cooper pairs" are extremely small with 
coherence lengths comparable with the size of the unit cell. Various mecha- 
nisms can lead to this local pairing HJ. In this paper we consider a model which 
can be relevant to the high-T c superconductivity because it contains not only 
such a local pairing but also an on-site electron-electron repulsion, interaction 
which may lead to the insulating phase of the cuprates. 

The Hamiltonian of the Penson-Kolb-Hubbard (PKH) model is [[J 

H = -tE (4+i,aCi,<j + H.c.) + U E n^n itl 

i,a i 

(1) 

+ WE(clclc l+lA c l+1 ^ + H.c.) 

i 

where we have used the standard notation for fermion operators : c\ a (c ia ) 
creates (destroys) an electron of spin o =t? 1 on the lattice site i and rii t(T = 
c\ a c icT . In the absence of the on-site interaction term U , the hamiltonian (1) 
reduces to the Penson-Kolb (PK) model || where the competition between 
the single- and pair-hopping of electrons can lead to interesting effects in ID; 
a spin-gap transition at half filling for W < has been the subject of some 
controversy [0, £|, Recently, some variantes of the PKH model were 
solved by Bethe ansatz method || [7j], where the single-particle hopping term 
is modified to include interaction effects: the jumping of an electron to an 
empty site differs from that corresponding to an occupied one |§. But the 
integrability of such model is possible only under some restrictions on the 
interaction parameters. Consequently, it remains interesting to have results 
depending on all parameters of the model and in a wide range. That is why 
we will consider both positive and negative values for U and W (t > 0) in 
Eq. (1) and arbitrary density; however, our results for the ground-state phase 
diagram are valid in a definite range of parameters determined below. 

We investigate the possible occurence of instabilities in the ground-state 
of the PKH model in the same manner || as it was done for the ID (t, U, X) 
model [|Hj : within the (zero-temperature) Green-function formalism in the 
Bloch representation, the instabilities are signaled by the poles of the vertex 
function T which obeys the Bethe-Salpeter equation. We solve this equation in 
the approximation when the irreducible vertex part is just the bare potential 
and the single-particle propagator has the 'free' expression. The imaginary 
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part of the poles (in the total frequency variable), interpreted as the inverse of 
the relaxation time to a new ground-state, gives us the regions in the param- 
eters space where the instabilities can occur; in the regions common to more 
instabilities we choose that phase with the shortest relaxation time. 



2 Bethe-Salpeter equation 

To understand what kinds of instabilities can occur in the system, one 
can investigate the generalized susceptibilities. It is assumed that we start 
from a phase where there is no order parameter, and study the density-density 
fluctuations. When the generalized susceptibility is singular, it is an indication 
that a spontaneous distortion or ordering can occur in the system. 

The general form of the susceptibility is 

X (k, u) = -*J dt e iu}t < T{V(k, t)V\k, 0)} > (2) 

where T is the usual chronological operator and T>(k,t) is a density operator. 
The brackets indicate an average in the ground-state (for zero-temperature 
case) or a statistical average (for finite-temperature case). If T>(k,t) is a 
charge- density operator, the corresponding susceptibily T CDW will test the 
charge density wave (CDW) instability; in a similar way can be defined the 
susceptibilities relevant to the other kinds of instabilities: T SDW for spin den- 
sity wave (SDW), T ss for singlet superconductivity (SS) and T TS for triplet 
superconductivity (TS). The correspondence between these quantities and the 
components of the vertex function can be found in Refs. |J. 

The generalized susceptibilities are two-particle Green functions and they 
obey the Bethe-Salpeter equation. In the simplest approximation, where only 
the bare interaction is considered for the irreducible vertex part and the one- 
particle propagator G is replaced by the free one G°, it can be written as: 

£T{k, k 1 ; K, Q) = —V(k, k'; K)+Y] V{k, k"; K)Q(k"\ K, Q)T(k", k'; K, O) (3) 

where Y can be any of the quantities T CDW , T SDW or T ss ; the TS case does 
not occur in this approximation because the interaction in the Hamiltonian (1) 
is only between electrons with opposite spin. K (Q) denotes the transfer mo- 
mentum (frequency) in the particle-hole (ph) channel and the total momentum 
(frequency) in the particle-particle (pp) channel, 

t = f - 1 CDW (A) 

4 ~ 1 1 SDW, SS [ ' 
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j +00 

— / dcu G°(k + K/2, u + SI/2) x 
2ir -00 



[G°(k - K/2, u - Sl/2) CDW, SDW 



(5) 



[G°(K/2 — k,Q/2 — u) 



SS 



where the addition or subtraction of the A;-vectors are defined modulo 2rr (the 
lattice constant is considered one). V in Eq. (3) comes from the interaction 
part of the Hamiltonian (1) in the Bloch representation and has the expresion 

v(k jj.m- 1 / U + 2Wcos(k + k>) CDW, SDW 
V{k,k,K)--<^ u + 2 Wcos(K) SS [b) 

where N denotes the number of sites in the chain. 

For the ph channel, Eq. (3) admits a solution of the form 



2ixN 



E T (k)X(K,n)E(k') 



with 



E(k) 




, x(K,n) 



X n 


X12 


X 




X22 


X 


X31 


X32 


x 



(7) 



(8) 



where E T means the transposed matrix of E. The unknown coefficients 
Xij(K,Q) are determined from the following algebraic system 



MX 



( U 


V 





-2W 




\ 


-2W J 



(9) 



where M is the 3x3 matrix 



M = 



f i-gU - Cl U 

-2 Cl W i-2c 2 W 
\ 2 Sl W 2pW 



SiU \ 
-2pW 
i + 2s 2 W J 



(10) 
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and 



iv q 



c n = 77 E (cos q) n G{q; K, ft) , n = 1, 2 

iv g 



Sn = 77 E (sin g) n £(g; if, ft) , n = 1, 2 

iv g 



P = 77 E sin q cos g£(g; ft) 

iv g 



(11] 



Since for the SS case k'; K) = V(iT), the solution of the Eq. (3) in the 
pp channel reads immediatly 



V(K) 



ss 



2nNl- V(K)g(K,Q) 



(12) 



An instability in the ground-state of the system occurs when the determi- 
nant D of the M matrix vanishes for the first time starting from the nonin- 
teracting case, indicating that the diverge and consequently T diverges. 
Following Refs. 0, we look for the T-poles of the form 



ft = E exc + iT 




CDW,SDW (K = 2k F ) 
SS (K = 0) 



(13) 



where E exc is the excitation energy to provide the system to undergo a phase 
transition; T is the inverse of the relaxation time of the unstable ground- 
state and can be also regarded as a 'temperature'; k F = ttti/2 is the Fermi 
momentum (n being the density of electrons) and e F = —2tcosk F . In this 
case, the determinant D has the form 



D ~ ji + p F X In 



^0 

T 



(14) 



which is valid for | T/ft |<C 1- This condition is similar to the BCS theory, 
where only the excitations of electrons around the Fermi level (with energies 
much less than the Debye energy) are taken into account. We will use the 
expression (|TJ]) to find the solutions of the equation D = for |T/ft | < 1 
and we expect the results to be reliable at least for not too big values of |A///| 
which plays the role of the coupling constant. 
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The parameters in Eq.(|H|) are given by 



2 w 2 1 -In \cosk F \ 

±1 + - (sin k F ) w ± — ± — —ri «w+ 

7T 7T 2 27T 2 siir fc F 




-f — In I cos k F \ + kp — In 2 I cos k F I ot 2 



SDW 
CDW 

ss 



(15) 



2 / In cos k F \ 2 2 

m + 2u> ± — sin K F H : uw ± — (sin k F ) w + 

TV \ SHI ftp / 7T 



A = < — 1-2 



+ 



kp 



tan fcp sin fcp 

. u + 2w 
where u = U/t, w = W j 1 ; 

fin = St sin 2 fcp < 



+ 2 In | cos k F \ uw z 



{coskpY 1 CDW, SDW 



SS 



SDW 
CDW 



SS 



(16) 



;i7) 



p F = (2tt sin fcp) (18) 

Pp/t being the density of states at the Fermi level. 

It follows that a transition to an ordered phase will occur at the critical 
'temperature' 

H \ 



| fin | exp 



(19) 



with /i/A < (so that |T c /fi | < 1). In order to get the phase diagram for the 
Hamiltonian (1), we determine at first the regions in the (w, u) -space where the 
quantity /i/A is negative (for each case: CDW, SDW or SS); when more than 
one instability can occur in a given region, we decide for the phase which is held 
first, i.e. with the shortest relaxation time (or equivalently, with the biggest 
critical 'temperature' T c ). We restrict our considerations to that region of the 
(u;, w)-space containing the origin u = w = and where A//x never becomes 
infinite. 
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3 Phase diagram of the PK model 



For U = the Hamiltonian (1) reduces to the PK model 0. By comparing 
the various critical 'temperatures' (for CDW, SDW and SS) we get the phase 
diagram plotted in Fig. 1: the investigated region is \w \ smk F < tt following 
from the condition |T c /f2 | < 1, as discussed above. For w > we get only 
a SDW. For w < there is a SS phase for densities less than a critical 
one n c > ^Arccos(e~ 2 ) ~ 0.914 and a CDW phase for n > n c ; the critical 
density n c tends to one as n — * 1. According to our calculations, at half 
filling the system is in a SDW state for w > and in a CDW phase for 
w < 0. Let us remark that at half filling the condition |T c /f2 | = 1 (when the 
effective coupling constant becomes infinite) determines the limits w = ±7r. 
It is interesting to note that close to our limit w = —ir, around the value 
w ~ —3.5, Sikkema and Affleck || found a phase separation transition; in our 
approach, the existence of such a transition can be in principle analysed by 
calculating the compressibility in the homogeneus phase [12| . 

Since at w = the electrons move freely and at large negative w the ground- 
state contains only doubly occupied and empty sites, it was argued that 
should be a 'pairing transition' at some negative w c . Exact diagonalizations on 
chains up to 12 sites [Q, |3|] have shown that at half filling this transition occurs 
around w c ~ —1.4; but conformal field methods [|J and renormalization group 
studies have shown that w c = 0. In our calculation the only transition 
which occurs at half filling (in the investigated range) is for w = 0. 



4 Phase diagrams for the PKH model 

In the limit of half filling (n — > 1) we found a SDW — CDW transition 
along the curve 

(20) 

in agreement with the prediction of Hui and Doniach M who found such a 
transition for u, \w\ 1 along the line u ~ — 2w. The obtained phase dia- 
grams at various densities are presented in Figs. 2-4. The first remark is the 
superconducting phase can not appear when u + 2w > 0; this is due to the form 
of the bare potential V(K), given by Eq. (6) in the pp channel, which becomes 
repulsive in that region. However, it follows from Figs. 2-4 that we do not 
get systematically a SS phase for V(K) < 0; some regions of SDW or (and) 
CDW phase still remain. The SS region decreases by increasing the density; 
it disappears for n — 1. Let us note that near half filling our phase diagram 
is qualitatively in agreement with that obtained by Hui and Doniach || who 



u 



-2w 



w 



l + (l-l/7r 2 )w 2 
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studied the same model with u > and w < using exact diagonalization for 
samples of up to 12 sites. For example, for rather big values of u and w they 
found a sequence SDW — > CDW — > SS in passing from u > —w to u < —w; 
this one can be also observed in our results from Fig. 2. 

5 Conclusions 

In this paper we have presented phase diagrams for the PKH model at 
arbitrary densities of electrons and for moderate values of the parameters W/t 
and U/t. Our mean-field- type approximation predicts results consistent with 
other works done at half filling. In the particular case of the PK model (U = 0) 
we have found for |W|/t < 7r/ sinfc F a phase diagram similar to that corre- 
sponding to the Hubbard model in the same approximation ; at half filling, 
the only transition which occurs is a SDW — CDW at W = 0. However, 
beyond the limits = ir/ sink F indicated by our approach, we expect 

a qualitative change in the ground-state of the PK model. To what extend 
this fact can be related or not with the phase separation transition found by 
Sikkema and Affleck at W/t ~ —3.5 near half filling, is a subject for further 
investigations. 
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